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Abstract 

Using loop equations, we compute the large deviation function of the maximum eigen- 
value to the right of the spectrum in the Gaussian (3 matrix ensembles, to all orders in N. 
We then give a physical derivation of the all order asymptotic expansion of the right tail 
Tracy-Widom (3 laws, for all /3 > 0, by studying the double scaling limit. 
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1 Introduction 



Extreme value statistics is an important issue in probability theory, which has a lot of 
applications in various fields, such as geology (earthquakes, etc), meteorology (extreme 
floods, etc) and hydrology, the question of establishing safety norms for construction 
works, growth problems in physics, finance, etc. Fischer and Tippett [FT28], completed 
by Gnedenko [Gne43], have shown that the distribution of the maximum of a large 
number of independent identically distributed random variables falls in one of the 
three universality classes known as WeibuU, Frechet and Gumbel, depending on the 
law of the initial variables (cf. also the book by Gumbel [Gum58]). However, the 
classification of limit laws for the maximum of a large number of correlated variables 
is a hard and still open problem. A possible strategy is to deal with particular models 
of correlated random variables, determine the limit law for the maximum, and try to 
show its universality within the model. 

The eigenvalues of random matrices provide a good example of such highly cor- 
related random variables. For the Gaussian ensembles, ie for N x N matrices with 
independent Gaussian entries, the distribution of the maximal eigenvalue is known to 
converge in the limit — )■ oo to the famous Tracy- Widom laws. Since the seminal 
work by Tracy and Widom [TW94, TW96], these laws have appeared in very different 
and surprising contexts, ranging from statistical physics and probability theory all the 
way to growth models [PSOO, TSIO] and biological sequence matching problems (for 
reviews see [AD99, TW02, Dei07, MajOT, ForlO, KKIO] and references thereof). 

For the usual Gaussian matrix ensembles GOE (Gaussian Orthogonal Ensemble), 
GUE (Gaussian Unitary Ensemble) and GSE (Gaussian Symplectic Ensemble), the 
limiting law of the maximal eigenvalue is given by the Tracy- Widom laws with Dyson 
index /3 = 1,2 or 4 (resp. for the GOE, GUE or GSE), which have been studied a 
lot. One can also consider more generally Gaussian beta ensembles for /3 > (G/3E). 
In this case also it has been shown that the maximal eigenvalue converges to a beta 
Tracy- Widom law TW^(s), but very few is known about these general laws. The goal 
of our paper is to study in details the asymptotics of these laws for s — oo (right tail). 

For that purpose we actually compute first the large deviations of the maximal 
eigenvalue in the G/3E, which describe its large fiuctuations on the right of the mean 
value (rare events). We expect that the large deviation matches smoothly the Tracy- 
Widom law, which describes the small typical fiuctuations around the mean value. By 
zooming close to the mean value in the expression of the right large deviation, we thus 
get the right asymptotics of the Tracy- Widom laws for /3 > 0. 

In the first section, we review some known results about the Tracy- Widom beta 
laws, for /3 = 1, 2, 4 and then for general /3 > 0. In the second section, we explain our 
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method and give our results about the right large deviation of the maximal eigenvalue 
in the G^E and the right tail of the Tracy- Widom beta laws. 

1.1 Usual Tracy- Widom laws 

In random matrix theory, the Gaussian ensembles are one of the first proposed and most 
studied ensembles. A Gaussian random matrix is a x matrix $, which belongs 
to some vector space Ejs, and whose entries are independent, identically distributed 
Gaussian variables. The distribution is proportional to d$ e~^"^'^^^'^^ with V^($) = 
and is invariant by conjugation under a group Gf^. (3 G {1, 2,4} is called the Dyson 
index: /3 = 1 defines the ensemble named GOE, for which Ei is the space of real 
symmetric matrices and Gi is the orthogonal group ; /3 = 2 defines the GUE, for which 
E2 is the space of hermitian matrices and G2 is the unitary group ; /3 = 4 defines the 
GSE, for which E4 is the space of quaternionic self-dual matrices and G4 the unitary 
symplectic group [Dys62, Meh04]. The distribution induced on eigenvalues is then 
given by: 



Because of the Vandermonde determinant A(Aj) = ni<j('^j ~ appearing to a power 
/3 in the distribution (1-1) for G^E, the eigenvalues are strongly correlated random 
variables. The distribution of the maximal eigenvalue when A^ is large is thus of special 
interest, and is very different from the possible limit distribution for the maximum of a 
large number of independent random variables, which are Gumbel, Frechet or WeibuU 
distributions. 

The fluctuations of the maximum eigenvalue in G/3E follow when A^ is large, upon 
rescaling, the celebrated Tracy- Widom laws. Thanks to orthogonal polynomials and 
integrable properties, these laws are very well understood when /3 = 1,2 or 4. Tracy 
and Widom have indeed shown more than 15 years ago the following theorem: 

Theorem 1.1 [TW94, TW96] In the G^E for (3 G {1,2,4}, converges m law 

to a random variable S distributed as: 



d/^A^(^i' . . . , Ajv) = W dAi e 



l<i<j<N 



(1-1) 



ns <s) = ^p{s) 



(1-2) 



Let q{s) be the Hastings- McLeod solution of Painleve II [HM80]: 





(1-3) 
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Define R{s) = q\t) dt = {q'{s)Y - sq\s) - q\s) and: 

E(.) = exp ( - - q{t) dt) , H(s) = exp ( - - R{t) dt) (1-4) 

Then: 

ci = l Fi(s) = E(s)H(s) (1-5) 

C2 = 1 F2(s) = H\s) (1-6) 

C4 = 2-^1' F4(.) = ^(E(.) + E-H.))H(.) (1-7) 

Furthermore, Deift and Gioev [DG07] have proved that, for polynomial V and when 
the large spectral density behaves as a squareroot at the edge of the spectrum (this 
is the generic behavior), the fluctuations of A^ax are still of order N'"^/^ and follow (up 
to rescaling) the Tracy- Widom laws. Actually, they establish universality at the edge 
of the spectrum for the correlation densities of n eigenvalues, for any n > 1. Those 
statistics can thus be expressed as in the G/3E, in terms of the Airy kernel [Meh04]. 

1.2 (5 Tracy- Widom laws 

The expressions in Thm. 1.1 make clear that the distribution of the rescaled fluctuations 
of Amax is very sensible to the value of /3, and one may wish to study the probability 
density in Eqn. 1-1 for arbitrary /3 > 0. Dumitriu and Edelman [DE02] have shown 
that Eqn. 1-1 with arbitrary /3 can be realized by eigenvalues of a random tridiagonal 
matrix, whose entries are independent but accurately chosen random variables. Given 
the lack of known integrability properties when /3 ^ {1,2,4}, this connection with 
random tridiagonal matrices and thus to a diffusion process has been the main source 
of results. 

Rider, Ramirez and Virag have shown the existence of a Tracy- Widom /3 law with 
the same rescaling: 

Theorem 1.2 [RRV06] In the G/3E for any /3 > 0, ^'^2/3^ converges in law towards a 
random variable S. The Tracy-Widom (3 law can be defined as its repartition function: 

E[S < s] = J\Np{s) (1-8) 

Notice that the definition of TW^ differs from the definition of F/3(s) when /3 = 4 by a 
rescaling of the variable s, i.e. TW4(s) = F4(2^/^s). Later, Virag and Bloemendal have 
provided a characterization of TW«: 
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Theorem 1.3 [BVIO] There exists a unique bounded solution to the problem: 

ldl + {s-W^)d^ + ds\M^^s)=0 

lim^„^_oo s) = 

liiiis^+oo limu,^+oo 0/3 (w, s) = 1 

and T\Nfs{s) = lim^^+oo (ppiw, s) . 

This can be reformulated in terms of a Schrodinger equation with time- dependent 
potential. Indeed, if we set: 



w, s ) exp 



4 V 3 



sw 



h{s) 



:i-io) 



1pl3{w,s) = 

for a choice of h{s) yet to determine, we find that ip/^^w, s) is solution to: 
with the time- dependent potential: 

However, it is not easy to perform analytic computations with TW^(s) starting from 
this characterization. One would prefer, for instance, a genuine differential equation 
characterizing T\Nf^{s) itself. 

For /3 = 2, it is possible to realize ip2{'w, s) as the solution of a 2 x 2 Lax system: 



P / 2 \2 , 

— (w - s) + 



l]w-h'(s) 



[1-11) 



[1-12) 



"^{w, s) 



[1-13) 



'ihiw,s)\ ( dyj'^{w,s) = L{w,s)'^{w,s) 

^2(^^,5) y' \ ds'^{w,s) = M{w,s)'^{w,s) 

The matrices L and M correspond (if w were rescaled by a factor ^) to the Lax pair 
of Flaschka and Newell [FN80]: 



Lfw, s) 



2 2 
—wq{s) 



q\s) 
- Q'is) 



-wq{s) 



+ q'{s) 
+ q\s) 



M{w,s) 



-f qis) 
qis) f 



(1-14) 

if we make the choice h'{s) = R{s) = {q'{s)) — sg^(s) — q^{s). The compatibility 
condition of the system 1-13 is the Painleve II equation for q, and for any s, there 
exists a unique solution subjected to the asymptotic condition: 

ip2{w, s) eo exp I - ^ - sw 

and the normalisation lims_j._|_oo lini^_^+oo V^2(wi, s) = 
function r(s) of this system can be represented as H(s)^ in terms of the Hastings- 
McLeod solution of Painleve II, up to a constant factor (see e.g. [FWOl]). This is 
indeed the representation of TW2(s). 

However, when /3 7^ 2, it is not possible to find a 2 x 2 Lax system leading to the 
Schrodinger equation 1-11 for the wave function ipi3{w, s). It is also not clear if TW/3(s) 
is a tau function for an appropriate integrable system. 



w — > —00 (1-15) 
1. It can be shown that the Tau 
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1.3 Right tail expansion 

In this subsection, we review the known results for the asymptotic behaviour of TW/3(s) 
when s — ?■ +00. The expansion when s — ?■ +00 is well known for /3 = 1,2,4. For a 
general /3 > 0, there are much less results. The goal of our paper is to show how to 
compute recursively all the terms of the large s expansion of TW^(s) for all /3 > 0. 

When /3 = 1,2 or 4, the asymptotic expansion of TW^(s) when s — > +00 can 
be derived in a straightforward way from its representation in terms of solutions of 
Painleve II (Theorem 1.1). We find that: 
4.3/2 



R{t) dt 



q{t) dt 



35 ,/2 3745 o 805805 289554265 ^ 

167rs3/2L 23-3 27.32 210-34 215.35 

31241084875 ^15/2 23604769513325 „9 _2l/2^ 
218 ■ 36 ^ . 38 ^ +0{s J 

41 „3/2 9241 _3 _ 5075225 „9/2 5153008945 .g 



e 



2v/is3/4 



24.3 '29-32 213.34 219.35 

1674966309205 3985569631633205 „ 21/2^ 



223 . 36 228 . 38 

Theorem 1.4 The following asymptotic expansion when s — )■ +00 holds: 



TWi(s) 



e 



41 ,,0 9241 , 5075225 5153008945 



40Fs3/4 



1 - TWafs) 



e 



24.3 '29.32 213.34 219.35 

1674966309205 3985569631633205 ^ 

223 .36 ^ 228 . 38 ^ ^K-^ 

3 /2 

35 3745 . 805805 ./o 289554265 . 

23 . 3 27 . 32 210 . 34 215 . 35 



167rs3/2 



1 -TW.fs) 



31241084875 ^^^/^ 23604769513325 „9 _2i/2n 

218 . 36 ^ ~^ 222 . 38 ^ +0{s ) 

3 /2 

143 41509 , 20443229 a/o 15418569025 



5127r s3 



3/2 , -^^-^"-^ -3 _ ^"^^'^^^■^ -9/2 , ^■^•^^"^u^w^u g 

24 . 3 ^ 29 . 32 213 . 34 219 . 35 

3330409204735 .r,,. 4908974519795465 „ 21 /2x 



223 . 36 228 . 38 

When /3 is arbitrary, taking advantage of the connection to a diffusion process in 
a time-dependent potential, Dumaz and Virag have established the first two terms of 
the asymptotic expansion when s — > +00: 

Theorem 1.5 [DVll] For /3 > 0, the following asymptotics hold when s — )■ +oo.- 

1 - TW,(.) = .-¥+o[ans)-/^l (1.16) 
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In this article, we study heuristically the asymptotic expansion of TW^(s) when s — )■ 
+00. Our result takes the form: 



Proposition 1.1 (Heuristics) 

r(f) 



1 -TW.Is 



3/3 2/33'^/^ 

s * e 3 



(4/3)5 27r 

where Rm are polynomials of degree atmost m 
the first few are: 



exp 



m>l 



3m 
2 



:i-i7) 



Ri{X) 

R2{X) 



24 

5 

64* 



+ 9X- 39) 
2 19X + 36) 



1ml, wii/i rational coefficients, and 



'IIX 



4608 



-1105X^ + 3240X^ - 23325X' 



(1-19) 

34938X - 41433) (1-20) 



In particular, the first two terms are in agreement with the rigorous result of Dumaz 
and Virag, see Eqn. 1-16. For /3 = 1, 2 and 4, the orders written above also match with 
those in Theorem 1.4. In section 2, we present a method to find the Rm recursively. 
For this purpose, one only needs to write down the solution of a recursive system of 
linear equations (and no inversion of linear operator is required), and then determine 
its behavior in some scaling limit. 



1.4 Left tail expansion 

In this subsection, we recall for completeness the situation for the asymptotic behaviour 
of TW^(s) when s — )■ —00. 

When /3 = 1,2 or 4, this asymptotic expansion can be obtained again from the 
known asymptotics of q{s), except that results about the total integral of q{s) are 
needed to obtain the constant prefactor. This last step was performed by Baik, Buck- 
ingham and diFranco [BBd08], and we have: 

Theorem 1.6 



TWi(s) 



2-11/48 gC'(-l)/2 



exp 



Inlsl 



-3/2 



3^2 16 24^2 



TW 



2{S 



JVJJs) 



2^/24 eC'{-i) exp 

2-37/48 gC'(-l)/2 



128 1152v^ 
|sp ln|s| 3|sh'^ 
"12" ~ ~ 



+ 



exp 
3|s 



64 
|s|3/2V2 



+ 



6 3 

731 s| -9/2 



+ 



512 9216v/2 



16 

63|s|-6 
8192 



512 
63|s|-6 
256 

ln|s| Is 

+ 



+ 



■0{\s\ 
-3/2 



48 V2 
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Let us mention that for /3 = 2, the value of the constant prefactor has been 

conjectured originally by Tracy and Widom [TW94] and first proved by Deift, Its and 
Krasovsky [DIK08]. Notice that only integer powers of |s| appear in the expansion 
of InH(s) when s — ?■ — oo. Moreover, E(s) is of order 6xp(— -^1^) when s — )■ — oo 
according to [BBd08, Corollary 1.3], so is exponentially small compared to the term 
E~^(s) in TW4(s). For this reason, the asymptotics to all orders of TW4(s) are the 
same as that of TWi(s) except^ for a sign change in front of each half-integer power of 

ki- 
ln a previous work [BEMNll], we have shown heuristically that, for all /3 > 0, 
TW^(s) admits a full asymptotic expansion when s — > — oo, of the form: 



where are polynomials of degree m + 2 in |, with integer coefficients. Lm is 
reciprocal up to a sign: 

L^^) = i-^r L^{^) (1-21) 

The Lm can be computed recursively (although we do not have a straightforward 
relation of the form = /m(-^m-i5 • • • , Li)), their coefficients are related to the sym- 
plectic invariants introduced by Chekhov and Eynard [CE06], computed for the curve 
of equation y"^ = x + ^ — 2. The first few are: 

= -54 + 193X - 193X2 + 54X=^ 
L2(X) = 2(-1190 + 5247X- 8042X2 + 5247X^-1190X^) 

L3(X) = -119862 + 608591X- 1199970X2 + 1199970X^- 608591X^ + 119862X^ 
L4(X) = 2(-3467274 + 19186065X - 44032530X2 + 56724246X^ - 44032530X^ 
+ 19186065X^ - 3467274X^) 

The constant prefactor is given by: 



2 2. v^e^'(°^?) (1-22) 



where ^(s; a) is a double zeta function, namely the analytic continuation of the series 
defined for Re s > 2: 

X{s;a)= V — --, x{s; a-^) = a' x{s; a) (1-23) 

^ — ^ m + am' r 

m,m'>0 ^ ' 



^The coefficients in tlie asymptotic expansion of f^is) = TW4(2^/'^s), which is the usual definition 
of the Tracy- Widom function for /3 = 4, also differ by powers of 2. 
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Moreover, the reciprocity of implies a duality valid to all orders in the asymptotic 
expansion when s — )■ — oo: 

TW,(s) = (^)^(^^)-^fw.[(|)n] (1-24) 

where the tilde indicates that we take the other sign for squareroots of |s| in the 
expansion. Those heuristic results are in agreement with all previously known results 
on the Tracy- Widom laws. In [BEMNll], we have also argued that the coefficients 
of this asymptotic expansion are universal, up to a non universal normalization of s. 
The universality class is defined by the blow-up of the spectral density in presence of 
hard wall close to the edge of the eigenvalue spectrum, which is generically the curve 
of equation = x + + 2. 

The present article completes this work by studying with similar methods the right 
tail asymptotics of TW^(s) (i.e. the limit s — )■ +00). We postpone an argument for 
universality of the right tail expansion to a future version of this article. 

1.5 Method 

In this subsection, we explain how we can recover the right asymptotics (ie limit s — 
00) of the Tracy- Widom laws TW^(s) from the right large deviations of the maximal 
eigenvalue of the G/3E. 

We first study the large deviations of Amax, i-e. the behavior when — )• +00 of the 
probability density function of Amax^ 

SnA^") = -^a/i]C,/3[Amax > «] (1-25) 

when a is independent of and is located to the right of the bulk spectrum. Deviations 
of Ajnax of order one from its mean value a* = lim7v-s.oo /^iv,/3['^max] are rare events 
when is large: the probability ^Ar,^(a) is actually exponentially small. Under some 
assumptions, we will show that Qn,i3{cI') admits a large asymptotic expansion of the 
form: 

g^^pia) = Cn,p exp ( N'-' (1-26) 

fc>0 

Then, if we substitute a = a* + N~'^^'^s in Q^^\a), ie if we zoom close to the mean value 
(approaching it from the right) looking at fiuctuations of order N~'^^^ from the mean 
value a* (instead of order one), we will find that: 

gf\a = a* + N-^/^s) ~ N'-^ ^/'(s) (1-27) 

This means that all terms in Eqn. 1-26 become of order 0(1), which is not surprising 
since ^Ar,^(a* -|- sN'"^^^) must be of order 1 and not exponentially small when A^ is 
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large. We indeed expect that yu)^^[Amax < a] = 1 — /u]^^[Amax > cl] for a = a* + sN~'^'^ 
converges to the Tracy- Widom law TW/3(s). Actually, we will see that a stronger result 
holds: 

VA;>3, Qf{s) = M,(^^{s"'Y->^ (1-28) 

where M is a Laurent polynomial in |, with rational coefficients. 

This suggests, and this is the non rigorous point in our derivation, to interprete 
exp (^X]fc>o^fc(f) st^^^'^)^ as the asymptotic expansion of: 

TW' (s) = lim iV-2/3 g^^p[a = a* + N-^/^s) (1-29) 

Notice that Eqn. 1-28 holds only for k > 3. The terms for k = 0,1,2 have to be 
handled separately, and the details of our computations show that they combine with 
the large behavior of the constant N'"^^'^ Cn,i3 to give a finite result, which reproduces 
Eqn. 1-16 as well as a constant prefactor which agrees with the values known for 
/3 = 1,2 and 4. 



2 Large deviations 

In this section, we study the right large deviations of the maximal eigenvalue Amax in 
the G(3E, i.e. the maximum of the Aj where the Aj are random variables distributed 

2 

according to the law in Eqn. 1-1 with a quadratic potential V{x) = We keep t > 
arbitrary in the calculations, so that the result can be read directly with the various 
conventions of normalization found in the literature. Our approach is based on the 
so-called Schwinger- Dyson equations for correlators of eigenvalues. 

2.1 Beta ensembles and Schwinger-Dyson equations 

We want to compute ^]^^^(a), the probability density function of Amax, with respect to 
the measure defined in Eqn. 1-1: 

GN,l3i'^) ~ ~f^a/^]v,^[Amax > o] 
-mv(a) r 

^N,I3 J]-oo,a]iv-i l<i<j<N-l 

= (2-1) 



where: 



NV(x) _ 21n(a-a:) -r 

f oo if X > a 
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and where for an integer M and a function U the partition function Z^j ^ is defined as: 

poo 

^M,, = n/ dA.e-^^(^') n l^^-^.l' (2-3) 

i=l l<i<j<M 

In other words, can be computed in terms of the partition function Zj^^i^j^ of 

a (5 ensemble with A^ — 1 eigenvalues in a modified potential Vjsi^a depending on and a. 



For any probability measure defined as: 

M 

i 



dAi^,^(Ai,...,AM) = ^ndA.e-^^(^») n l^^-^.l'' (2-4) 

i=l l<i<j<M 



we can derive Schwinger-Dyson equations, which give relations between some expecta- 
tion values with respect to /i^/^j, by using integration by parts. To write them down, 
it is convenient to introduce the correlators: 



n M 



iy„(a:i, . . . , x„) = ftli,^ [(HE ^^v) J (2-5) 

j=iij=i ^ 



where c stands for cumulant and fJ'^.f i^l- ■ ■ ] denotes the mean value of ■ ■ ■ with respect 
to /i^^^. The Wn are related to n-point correlations of density of the Aj's. For instance, 
Wi is related to the spectral density 

i=l 

for X G M, in the sense of distributions. Let us also define: 



i=l * j=2ij=l ■> 



In those notations, it is implicit that Wn and P„ depend on M, /3 and U . We also use 
the notation x/ = (xj)jg/ when J is a set. The Schwinger-Dyson equations (also called 
loop equations) for the beta ensemble have been established many times [CE06, BGll]. 
We give them here without proof: 

Theorem 2.1 For any x G M \ f/^^({+oo}).- 

W2{x,x) + {Wi{x)f |)(5^.Tyi)(x) -MU'{x)Wi{x) + M Pi{x) = (2-8) 
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And for n > 2, if I denotes the set {2, . . . , n}, and z, 0:2, . . . , x„ G M \ f/ ""^(l+oo}) are 
n points distinct two by two: 

+ ^ VF|j|+i(x,xj)VF„_|j|(x,X/\j) + (l - ^){d^Wn){x,xi) (2-9) 
jci \ p/ 

-MC/-wiy„(...,) + MP„(...,) + |^3..( ""-'"'"'^w'-'^'-<"'' ) = 



We stress that those relations are exact for any finite M . Unless explicitly stated, we 
shall apply these relations to the measure /iyv^i^^ only, i.e. we only use correlators 
with respect to the modified measure [ij^^^^^. 



If we know the correlators, we may come back to the derivative of the partition 
Lction. 
we have: 



function. Since the integration weight in yu)^^'^ o vanishes when one of the Aj equals a. 



dMZl''il^ = m^{a) (2-10) 

We then need to integrate the above expression to get In Z^j^^^^ ^ and then the probability 
density function of Amax, ie Q^^^ia) (cf Eqn. (2-1)). To obtain the integration constant, 
let us remark that for any fixed A^: 



zl^il ~ (2-11: 



NV 



Accordingly we get an explicit expression for the probability density function of Amax, 
ie pjySi) as a function of the first correlator W\. 



iV/3 



Q\,{a) = ^-^'^ a^^-)/^ exp - /3 / da'(l^,(a') - ^) (2-12) 

So, if we want to compute ^]^^^(a), our first goal is to compute Wi{a), and for this task 
we can use Schwinger-Dyson equations. 

2.2 expansion and large deviations 

Under reasonable assumptions on V (which carry on to VN,a), one of the author and 
Guionnet showed in [BGll] that the correlators with respect to fi^l\ p admit an asymp- 
totic expansion when — t- +00, of the form: 

PF„(xi, . . . , a;„) = 5^ Ar2-"-^ W(a;i, ...,x^) (2-13) 

fc>0 
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when points (not necessarily distinct) in a neighborhood of oo which 

does not depend on A^. The important feature in Eqn. 2-13 is that Wn is a 0(A^^~"). 
Let us plug this expansion in the Schwinger Dyson equation. But notice first that: 

(A^ - 1) VN,a{,x) = N V{x) - 2 ln(a - x) (2-14) 

allows us to decompose: 

{N-l)Pnix;x2,...,Xn) = NQn{x;x2,...,Xn)-'^^^ (2-15) 



where: 



i=l * j=2 ij=l J 



The large expansion for the correlators also implies an asymptotic expansion for Qn, 
of the form: 

Qn{x- X2,...,Xn) = Y, Qn'(^; ^2, • • • , (2-17) 

fc>0 

At leading order for Wi, we find: 

(lyfl (x))' - r'(x) Vrfl + Q^?{x) = (2-18) 

which implies that W^j^^ does not depend on /3 and a: it coincides with the leading order 
of the correlator with respect to /i]^^^25 i-^- of the measure for a random hermitian 
matrix with potential V. For the subleading orders, we find: 



{V'{x) - 2Wf\x))wP{x) - Qf\x 
Wt'\x,x) + J2wl'\x)wt 



k'=l 



and for n > 2: 



l-|)awf-')W.2 '^^"W-'^-""M (2-19) 



{V'{x) - 2Wl°\x))Wjt\x, xj) - g|fl(x; xj) 

Wl'-,'\x,X,Xj) + Yl ^|J|+l(^'^^) Wn-\J\i^^^I\j) 



JCI 0<k'<k 
(J,fc')7^(0,O),(7,fc') 



pJ X — a 

/3 ^ ' V X - Xi / 



is/ 
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In fact, Qn{x',xi) = {O ■ Wn){x,xi) is found by action of a linear operator O on 
Wn- Provided one knows how to invert [V — 2H/"|°')ld + O, the relations above form a 
triangular system which can be solved recursively to find Wn{xi, . . . , Xn)- The recursion 
is on the integer n + k. 

2.3 The Gaussian case 

2 

When V is a. quadratic potential, i.e. V{x) = |j, it is even simpler, because all the 
quantities Qn are known a priori: 

- 1 

Qi{x) = 6n^i — - — , Q„(x;x2, . . . = for n > 2 (2-21) 

6 



Moreover, when a > 2\/t ie a > a* (for V{x) = %r the mean value of A 



IS 



a* = 2\/t when — oo), the effective potential V^f^ satisfies the assumptions: 

o V)^f^ : ] — oo, a] — 7- M is continuous, and admits a 1/{N — 1) expansion. 

o its large (A^ — 1) limit is V'^^^{x) = and it is known that the large (A^ — 1) 
spectral density with respect to /i^'^^ ^ has a connected support It = [—2\^, 2\^], 
and behaves as a squareroot ±2y/i. 

o When a > 2\/i, V]\f^a is analytical in a neighborhood of It. 

According to [BGll, Proposition 1.1], this ensures the existence of a 1/(A^ — 1) expan- 
sion for the partition function and the correlators, which can be repackaged into a.l/N 
expansion as stated above in Eqn. 2-13. 

Eqn. 2-21 imply that, for the Gaussian potential, Eqns. 2-19 and 2-20 really form 
a recursion on n + k without any linear operator to be inverted. The leading order of 
W^i (cf Eqn. 2-18): 

W!>\.) = (2-22) 
is associated to a large A^ spectral density p[x) given by (cf Eqn. (2-6)): 

- ^ (2-23) 

We recover without surprise the Wigner semi-circle law, which is valid for all values 
of /3 > [Wig58, Joh98]. Let us write Y{x) = V'{x) - 2Wl^\x) = \Vx^ - At. To 
compute the subleading orders, it is convenient to introduce a uniformization variable: 
z, which maps xeC \ [-2v^, 2^] to C \ D(0, 1): 

. = yi(. + -) « (2-24) 
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and to work with modified correlators: 

a;lf](zi, ...,Zn) = x'(zi) ■ ■ ■ Wn{x{zi), x{zn)) (2-25) 

These quantities are now rational functions of Zi,...,Zn- Let us call a such that 
x{a) = a and |a| > 1. One can prove by recursion that the oon^ have poles only at 
Zi = 1,-1, a~^, and Zi = zj^ for i ^ j. 

From Eqn. 2-22, we have ujf\z) = {z~^ — z~^), and we find at the first step of the 
recursion, cf. Eqn. 2-19: 



Then, at the second step, cf. Eqn. 2-20: 

^l'(-) = ;^2.: in2 (2-28) 



(a2 _ 1)2(2 -a-i) 

5(f)'-9f + 5 (f)%^(-f + l) -(D^I + ^^^S^ 

16(z-l)4 ^ 16{z-l)3 32(z-l)2 

c,^'2^2 ,q2 r /2\2 I 9a-7 ( 2 , i ^ /- 2 'i 2 2 , -7a^ + 18a-7 



16(z + l)4 16(z + l)3 32(z + l)2 

^ iG'^)[ (..-i)i2-i)- - (..+i)3i2+i)^ 

2z2 ^|(l + 3z|) 



(Z1Z2-1)3(Z2-1)2 1)2(^1-1) 

1 2 1+a 1 2 l-g 

2 1-a 2 ^ 1+a 



(Zi - 1)2(Z2 - 1)2 (Z1 + 1)2(Z2 + 1)2 

[01. N 8 {I + ZiZ2 + Z2Z3 + Z3Zi){zi + Z2 + Z3 + Z1Z2Z3) /o on^ 

-3(-l--2,Z3) = ^ (,2_i)2(,2_i)2(,2„i)2 (2-30) 

and so on. Despite its non-symmetric expression, 00^2^^!^ ^2) is actually symmetric in 
zi, Z2, and in general Un\zi, . . . , Zn) must be symmetric functions of all their variables. 

Actually, we remark that the auxiliary quantity Pn{x]Xj), appearing in the 
Schwinger-Dyson equations 2-19 and 2-20, is a priori explicitly known in all the clas- 
sical ensemble, i.e. Gaussian, Laguerre, and Jacobi ensembles. This simplification 
makes the approach based on Schwinger-Dyson equations attractive to study large 
expansions in those ensembles. 

2.4 Large deviations of the maximal eigenvalue in the Gaus- 
sian case 

Using Eqn. 2-12 and the expressions obtained recursively above for uj^^\z), we can then 
find after integration the desired asymptotic expansion (Eqn. 1-26) of the probability 
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density function of Amax, ie ^^'/f (a)- The prefactor in Eqn. 2-12 can be computed 
exactly using Selberg's integrals: 



2 , 3-1-/3" 



X exp ( - — — ) cxp 



t^-^a^''-''^ (2-31 

fc>0 



^^/? A r o V- ..i-fe r ,W ^./^ _ ('^M-4,i)(«^^-l) 

a'(a'2 -f 1) 



where we recall a{a) = "' 2'^ ie a = a/^ (a -|- ^) . More explicitly, using Stirling's 
formula: 

Proposition 2.1 G^'^^i^) — ~dafi'^^^[\ma.x > a], for a > 2y/i, has a large N expan- 
sion of the form: 

-4 (l\'^-^r(§. 



In 



27r I Ut ^ \ 2Vi 



2V 2;V^/V2 /\ 2V^ 

m /.oo 



where -Bm+i ore t/ie Bernoulli numbers (and they vanish when m + 1 is odd). 

This large asymptotic expansion describes the probability of a large deviation of 
Amax to the right of the edge of the spectrum. The leading order is in agreement with 
the result of Majumdar and Vergassola who computed it using a Coulomb gas method 
[MV09], and the first correction has been found recently for /3 = 2 by Majumdar and 
one of the authors using a method of orthogonal polynomials on a semi-infinite interval 
[NMll]. We give below the three first decaying corrections: 
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r , [21 , -5| +27|-39 -3|f + 19|-33 |-4 

/a(a) 6(a2_l)3 4(^2-1)2 2(^2 _ 1) 

r , ' [3], -10(|)^ + 73(f)^-19lf + 180 

-25(|)^ + 187(|)^ - 507 1 + 501 -80(|)^ + 627(|)^ - 1807 1 + 1926 
^ 2(a2 - 1)5 ^ 8(a2_ 1)4 

-15(f)^ + 133(f)^^438f + 539 3(f)^^20f + 38 

6(a2 _ 1)3 + 4(c,2 _ 1)2 y --^ i 



7a(a) 



, [4], -1105(f )^ + 9720(f )^ - 34557(f )^ + 59238 f - 41433 
= 18(a2 - If 



985(|)* + 8724(|)^ - 31389(|)^ + 54786 1 - 39273 
4(a2 - 1)8 

767(|)^ + 6871(|)^ - 25157(|)^ + 45003 1 - 33321 
2(a2 _ 1)7 



-3443(|)* + 31476(|)^ - 118455(|)^ + 219640 1 - 170091 
12(a2 _ 1)6 

-1014(1)" + 9660(|)^ - 38180(|)^ + 75015 1 - 62150 
10(a2_ 1)5 

-105(|)* + 1120(|)^ - 4953(|)^ + 10902 1 - 10153 
8(a2 - 1)4 

15(1)' -128(1)' +412 1 -506 



2.5 Scaling regime and Tracy-Widom laws 

As explained in § 1.5, we expect that if we plug naively the scaling a = \/t{2 + N~'^/'^ s) 
in the large deviation of the maximal eigenvalues (Proposition 2.1) (where it is not 
valid as the large deviations describe fluctuations of order one from the mean value 
a*), we should recover the right tail of the Tracy-Widom laws (that describe small 
fluctuations of order N~'^^^ close to a*). 

We first need to find the behavior of ujf^\a) at leading order when a — )■ 1 (i.e. 
a — )• a = 2Vt). The cases k = and k = 1 can be computed from the expressions in 
§2.4, and in general we have: 

Lemma 2.1 For any k > 2, there exists a polynomial Rk-i of degree k, with integer 
coefficients, and a positive integer pk_i, such that, when a — t- 1.- 

^S" («) = o Tfr'^^l. 2 + O f 7 \^) (2-36) 
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The proof of this technical lemma is given in § 2.6. Since ooi can be computed recur- 
sively, so do the polynomials Rm- The first few can be read from § 2.4: 

pi = 3 i?i(X) = -5X2 + 27X-39 

P2 = 5 R2{X) = 3{-10X^ + 73X^ - 191X + 180) 

P3 = 9 RsiX) = -1105X^ + 9720X2 - 34557X2 + 59238X - 41433 

If we write a = \/i{2 + X"^/^ s), we have (a — 1) ~ ^/s when X is large, and combining 
the former lemma in Proposition 2.1, we find: 

Proposition 2.2 (Heuristics) The probability density of the Tracy-Widom law has the 
following asymptotic expansion when s — t- +oo.' 

TW'(s) = hm ytX-2/3^^'l'^(a = v^(2 + X-2/3s)) 

3 V2 4/ ^ 2P--3-m 



(4/3)^ TT 



rn>\ 



The polynomials Rm in the result announced in Proposition 1.1 are related to Rm by 
the triangular system: 

^'"(f) _ p ^-^ 2 ( - 1 + Ei>i^O! ( 3x j]^^^™, -p-r r .2..„ 
2."^.3-m-^™y- 21- ^ ^"4^ " llK^-l)^-i(^)J 

' m,>l (i>l) -^^^^l * i>\ ' 

with the convention {i — l)Ri\.__^ = 1. We have checked that this proposition gives the 
correct results for the orders written above when /3 = 1, 2 and 4. 

2.6 Proof of Lemma 2.1 

After successive partial fraction expansions, Un\zi, . . . , Zn) can always be decomposed 
in a unique way as a sum: 

[fc] / N _ \ ^ Mr{f3; a;zi,..., Zn) 

(2-37) 

where M^{[i\ a; Zi, . . . , Zn) is a rational function of its arguments which does not vanish 
when one of the monomial in the denominator vanish. It may have poles when Zi —1 
or 0, when a — )■ 1, but we are not interested in those poles. We define the degree of each 

term as 5{r) = XlILiiil'") + ^^l^") + ^^l^") + Z]i<i<i<n "^ijl'^)' ^^"^ define the degree 
of ' as dn,k = maxr S{r). The interest in this notion of degree is that Un\a, . . . ,a) 
must be a O ((a — 1 )""'".'=) when a — )■ 1, and in particular: 



uj\ [a 



G 0((a- l)"^^-^)) (2-38) 
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We also want to trace back the dependence of the terms of higher degree in /3. 
Eqns. 2-19 and 2-20 imply that wlf' is a polynomial in the variable |. Let bn,k, the 
maximum power of | appearing Mr{(3, a; Zi, . . . , Zn) when the degree 5{r) is maximal, 
i.e. is equal to dn,k- 

We now determine bounds on dn,k and hn^ki by recursion on n + k. Lemma 2.1 is 
then an application of the lemma below to n = 1. 

Lemma 2.2 For all {n,k) ^ (1,0).- 

dn,k < 4n + 3A; - 6, hn,k <n + k-l (2-39) 

Proof. When n + /c = 2, we already know the answer from Eqn. 2-26 and 2-27, namely: 



= 1 
hiA = 1 



0?2,o — 2 
^2,0 = 1 



(2-40) 
(2-41) 



and it satisfies the desired property. Let p > 2 and assume the lemma whenever 
n + k < p. When n + k = p + 1, Eqn. 2-19 or Eqn. 2-20 computes cun ' in terms of u^^, ^ 
with n' + k' < p. This relation can be organized as follows: 



uj}!:\z,zi) 



Y{x{z))x'{z) 



A^^_^_^\z; zj) + -Bi+i(2; zj] 



(2-42) 



where: 



AJ+i^1(2;2;j) = u^^^^'\z,z,zi) 



B 



[k\ 

n+l 



D^:Kz;zj) 



E^^l-^{z,Zi;zj) 



[k+2] . 



[k] 



[k-k']. 



JCI 0<k'<k 
(J,fc')^(0,O),(/,fc') 



x'{z) dz 



UJri 



x'{z) 



{x'iz)Yz ^'(^^ 



x'(a) 



z — a){z — a ^) 



(2-43) 
(2-44) 

(2-45) 
(2-46) 
(2-47) 



The very definition of the degree also allow to follow the degree of each of the terms: 

o The degree of A^^^^^ is atmost dn+i,k-2, which is by assumption smaller than 
An + Sk- 8. 
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o The degree of is bounded by the maximum of {dn'+i,k' + dn-n',k-k') for 

< n' < n, < k' < k and {n',k') ^ (0,0), {n,k). But since dn'.k' is bounded 
by an affine function of n' and k' for n' + k' < n + k hj assumption, all these 
numbers are bounded by the same quantity, namely {4{n' + 1) + 3k' — 6) + (4(n — 
n') + 3{k - k') - 6) = 4n + 3fc - 8. 

o The degree of Cn is dn,k-i + 1 < 4n + 3/c — 8. 

o A term in wjf"^' containing a factor {z — yield a term containing a factor 

[z — in Dn Thus, the degree of Dn ^'^ is bounded by dn,k-i + 1 < 

4n + 3A; - 8. 

o A term in cujfl^ containing a factor {zi — !)"•' yield in the ratio E^^^{z, zf, zj) 
several terms, with possible factors (z — l)^(2j — l)~^^~^^\zzi — {z — l^izi — 
l)-^^+^\zZi - 1), {z - l)2-0'+2)(22. - 1)^2 _ i)2-(i+3)(^^. _ lyi^ Tj^^g^ ^Yie 

degree of -D„_i,fe is bounded by dn^i^k + 2 = 4n + 3/;; — 8. 

Besides, the quantity Y{x{z)) x'{z) = —z^'^{z—l)'^{z+l)'^ has degree 2. So, we conclude 
that for n + k = p + 1, we have dn,k < (4n + 3/c — 8) + 2 = 4n + 3/c — 6 which is the 
desired bound. 

Let us come to the dependence in f3. We know that Y{x{z))x'{z) does not depend 
on /3. Eqn. 2-42 implies that 6„ ^ is bounded by the maximum among the numbers 
bn+i,k-2, K'+i,k' +bn-n',k-k' ioT < u' < u, < k' < k and {n',k') ^ (0, 0), (ra. A;), 
(6„,fc_i-|-l) and And the recursion hypothesis implies that all those numbers 

are bounded by n + /c — 1, hence the bound for 6„ fc. □ 

Eqn. 2-42 implies that a;lf' (2:1, . . . , z^) is a rational fraction in zi, . . . , 2;^, and also of 
the variable a, with rational coefficients. We would like to know which denominators 
do we get, especially in front of the leading term when a — )■ 1. Notice that Eqn. 2-42 
only involves integer coefficients, and the only poles in the partial fraction expansion 
of M,(/3, a; ^1, . . . , 2„) (cf. Eqn. 2-37) are of the form [a + l)--''^('^) VCi=x{'^i + l)"''-^'^- 
So, we can only get a power of 2 in the denominator, coming from expansion of [a -|- 
l)-^'(^-)' X\.{z + l)-*^^('^) when Zi a and a — ?■ 1. This justifies the factorization of a 
power of 2 (which we do not try to compute here) in Lemma 2.1, and the fact that Rm 
is a polynomial with integer coefficients. 

Addendum 

While we were completing this work, P. Forrester [Forll] has studied the two first lead- 
ing terms and the constant prefactor in the large deviation function, and its heuristic 
matching with the right tail of Tracy- Widom ji laws, with similar results. His method 
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relies on the functional equation satisfied by the spectral density and its first correc- 
tion. This corresponds in this article to the loop equation 2-8 up to 0(1), in which case 
the term W2{x,x) can be forgotten and the equation only involve the spectral density 
itself. [Forll] also obtain the two first orders and the constant prefactor of the large 
deviation function in the Laguerre ensemble. As we pointed out in § 2.3, the method 
of loop equations could be applied to obtain recursively, without any inversion of linear 
operator, the full asymptotic expansion of the large deviation function. 
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